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Abstract 



A theoretical model for magnetooptic imaging is presented. The model gives 
a detailed description of the magnetooptic indicator and the optical imaging 
system, where the intensity at the detector plane is described by a set of 
diffraction integrals. To demonstrate the applicability of the model, some 
issues concerning imaging of single vortices in superconductors are studied. 
It is found that the useful signal is strongly dependent on the penetration 
depth of the superconductor, and also the sensitivity and thickness of the 
indicator. 
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I. INTRODUCTION 



For centuries humans have tried to 'see' magnetic fields using various devices. Iron par- 
ticles, magnetic fiuids, inductive coils, Hall probes and magnetic bacterias are only a few of 
the approaches which have been adopted. During the past 40 years the research on mag- 
netooptic visualization has been highly motivated by the need to 'see' magnetic fields from 
superconductors. To that end, magnetooptic visualization using thin film indicators has 
been particularly useful. In 1957 P.B. Alers visualized for the first time the magnetic field 
distribution in a superconductor using a mixture of cerous nitrate-glycerol 0. A highly 
successful development followed in 1968 when H. Kirchner noticed that chalcogenide-films 
could be used to visualize fields with improved sensitivity 0. In the following 20 years 
magneto-optic studies were published periodically, mainly looking at fiux profiles in super- 
conductors 0,^. However, there were no major advances until 1991, when it was discovered 
that bismuth-substituted iron garnet films with in-plane magnetization are excellent in- 
dicators with almost no domain activity 1^-0. This discovery triggered a large number of 
quantitative investigations of the fiux behaviour in high-temperature superconductors |]8|-[T0 



Two of the most recent advances in this field are ultrafast magnetooptical studies (see e.g 



ref. [^) and imaging of single vortices |T2 



Although magnetooptical indicators have mainly been used to image magnetic fields 
from superconductors, a large number of other possible application areas exists. Thus, it 
has been realised that such indicators are useful for imaging magnetic defects, magnetic 



phase transitions, microelectronic circuits and magnetic storage media |T3- 

Since the field of magnetooptic visualization has matured for more than 40 years, one 
may expect that most parts of the image formation process is well understood. This is 
not the case. In fact, only recently a model for magnetooptic imaging has been published 
||18|| . However, this model is mostly concerned with the magnetization redistribution in 
the indicator. The purpose of the current paper is to develop a more detailed model for the 
magnetooptic imaging system (see Fig. including the indicator, and to use this to examine 
the possibility of imaging weak magnetic fields, e.g. single vortices in superconductors. 

II. MAGNETIZATION IN THE INDICATOR 

The basic sensor element for detection of magnetic fields is the magnetooptic indicator. 



It is therefore first necessary to know how the indicator respond to an external field [jT8[ . Let 
M be the magnetization vector in the indicator, and Mg be the saturation magnetization. 
Then we have from Fig. ^ 

M = Ms[cosPcos(f)M,cosPsin(f)M, sin(3] = [MrCos(j)M, MrSin(f)M, Mz] . (1) 
The external field is given by 

H = [HrCOSCpH, HrSin(j)H, Hz] , (2) 



where Mr = + and Hr = yH^ + H^- We assume that the indicator consists of a 

single domain, and that its free energy can be expressed as a sum of the uniaxial anisotropy 
energy (Eu), the Zeeman energy [Ez) and the demagnetizing energy (Efi) 
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E^Eu + E, + Ed . (3) 

The uniaxial anisotropy energy can be written as, 

E^ = K^cos'fi , (4) 

where is the uniaxial anisotropy constant. The Zeeman energy is given as 

Ez — —HqM-H — —^QMs[Hj.cos(5cos{(j)M — (pn) + HzSin(5] . (5) 

Since we neglect the cubic anisotropy of the indicator the in-plane magnetization direction 
must be the same as that of the external field (0m = 0if ) ■ The demagnetizing energy, which 
tends to keep the magnetization in the plane of the film, can be given as 

Ed = ^ sm (3 . (6) 
Minimizing E with respect to P results in the following relationship 

HaSinf3cosf3 + HrSin(3 — H^cosP — , (7) 
where the so called anisotropy field is given by 

2 

Ha^M, ^ . (8) 

If we assume that Hr — 0, then the resulting expression for the magnetization is 



\\i-[jf) ,H,<Ha , (9) 



and 



M, = M,^ ,H,<Ha . (10) 

-"a 

If /3 is very small (corresponding to very small H^, or large Ha or Hr), we may put sin/? ~ /? 
and cos/3 ~ 1, and the magnetization becomes 

Mr^M, ,H,^Ha , (11) 

and 

M, ^ M,— ^ ,H, « Ha , (12) 

Thus we see that an in-plane component Hr will only reduce the z-component of the 
magetization. Therefore, it is necessary to keep Hr and Ha small if large sensitivity to Hz 
is needed. 
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III. THE MAGNETOOPTIC IMAGING SYSTEM 



The basic imaging system is shown in Fig.|I]. In principle, it is a conventional polarizing 
microscope, although many different designs exists |T^. Basically, the z-component (or any 
other component) of the magnetization is converted into a small rotation of the plane of 
polarization, which can then be detected by an array detector (e.g. CCD). The signal on 
the detector can be found by an analysis based on Jones matrices together with the methods 
used in refs. [^,^. For a general polarization distribution incident on the objective we write 



Eq 



a(^,0) 
6(^,0) 



where 6 is the incident angle and (j) is the azimuthal angle. In this paper we assume that 
the incident polarization after the polarizer consists of a x-component only (a=l and b=0). 
The effect of the objective lens is to change the incident polarization into s and p-polarized 
light with the following transformation 

COS0 sinq 
—sincj) cosq 

Next, the rays passes through the isotropic substrate on which the indicator is deposited 

tf 

When the rays are reflected from the magnetooptic layer, some of the s-polarization is 
transformed into p-polarization (and vice versa). In order to assure maximum reflection, 
one may evaporate a thin mirror onto the indicator. The reflection matrix from the indicator 
and the mirror can be expressed as 



R 



^sp ^ss 



where r^p and r, 



are the usual Fresnel coefficients. In the first order approximation (cor- 
responding to weak magnetooptic effects), we find the following off-axis components for 
arbitrary magnetization vectors and angles of incidence 



ps 



TxriiQ 
An 



,p 
''01 ''10 



D 



'a{z)M^ + b{z)Mytan9i 



dz 



(13) 



r., = -'-^t'„A X M (^f^MMl^MM^^ i, , (14) 

where ni is the refractive index of indicator, Q is the magneto-optical material constant, 
Aq is the wavelength in vacuum, tlf is the s or p-component of the transmission coefficient 
through the interface from medium i to j (similar for r^^j^), 9i is the propagation angle in the 
indicator, f{z) = expl—iiTTnizcosOi/Xo], and 
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[l + rlJiD-z)][l + r{J{D-z)] 
[l-rfor'ufiDm-r'ior'ufm 



(15) 



b{z) 



[l + riJiD-z)][l~rlJiD-z)] 



(16) 



according to the recipy of Hubert and Traeger PT|. The integration is carried out over the 
thickness (D) of the indicator. Here we have neglected the transverse Kerr-effect, since this 
is rather small in most indicators. 

For strong magnetooptical effects, one should use a similar formalism, but take into 
account multiple scattering. This complicates the problem considerably, and will not be 
pursued here. 

On its way back to the array detector, the light will again pass through the substrate 



- out 





tr* 



Furthermore, it is collected by the objective lens, and in total we can write the electric field 
before the analyzer as 



E — L ^ToutRTinLEQ . 



(17) 



After evaluating the equation above we observe that also the Fresnel coefficients between 
the substrate and air influences the polarization state of the outgoing beam. In many cases 
it is possible to deposit a dielectric coating which minimizes the depolarization due to the 
substrate. Here we assume the most ideal case, in which all transmission coefficients become 
unity. In this case we may write 



and 



E^ 



E„ 



1 1 

r, Vpp + "^ss) + „ [T"pp 



rss)cos2(p - ^(r^p + rps)sin2(p 



2 i^sp fps) 



+ 2^'^PP - rss)sin2(t) + -(r^p + Tpjcos^ 



(18) 



(19) 



These equations can be used to calculate the expected polarization pattern at the exit pupil 
seen in the conoscopic image, which is often used to judge the quality of the imaging system. 
To find the electric field at the detector plane, (r',0c), we assume that the imaging system 
obeys the sine condition, and use the method of refs. p!9| , |20[] to obtain 



E^D = Iq+ /a cos20c + Ilsin2(j)c , 



(20) 



and 



E^ 



yD 



(21) 



where 
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^0 — {rpp + rss)sin9cos9Jo{kr'sin9)exp{2ikzcos9)d9 , (22) 
Jo 

pa 

^0 ~ (''"sp — rps)sin9cos9jQ{kr'sin9)exp{2ikzcos9)d9 , (23) 
Jo 

-^2 = / (rss — rpp)sin9cos9J2{kr'sin9)exp{2ikzcos9)d9 , (24) 

^ 

ra 

-^2 = / {^sp + rps)sin6cos9J2{kr'sin9)exp{2ikzcos9)d9 . (25) 
^ 

Here J„ is the Bessel function of the first kind (of order n), and a constant in front of 
the integrals is ommited. Note that the integral is taken over the angular aperture of the 
objective lens, which has a numerical aperture NA = sina. In deriving these formulas, 
we have assumed that the objective is corrected for all spherical aberrations introduced by 
the film and the substrate. Here the factor cos9 occurs since we have chosen aplanatic 
apodization at the exit pupil {\/cos9, squared for two passes through the objective). The 
phase factor exp{i2kzcos9) is included to account for the defocus of the light beam. Here 
we will assume that z=0, and therefore that the light beam is at focus on the mirror-plane. 
When an analyzer is placed in front of the detector, the intensity at the detector plane must 
be written as 

in = \E:rDC0S'y + EyDsiwyl^ , (26) 

where 7 is the analyzer's angle from the x-axis. So far, we have not considered diffraction due 
to the localized magnetization distribution in the indicator. However, this can be included 
by using the procedure of Kambersky et al. p2[|. In that paper a procedure for calculating 
the diffracted Kerr or Faraday amplitude is presented. The result must be weighted by the 
polarization vector of the lens system, and summed up at the exit pupil. This is a large 
numerical task, at least for a reflection-type microscope, and will not be pursued here. That 
is, we will neglect diffraction from magnetization gradients in the indicator. 



IV. IMAGING OF SINGLE VORTICES IN A SUPERCONDUCTOR 

Only recently magnetooptic imaging of single vortices was achieved experimentally [jl2 . 



In this section we study magnetooptic imaging of single vortices using approximative so- 
lutions of the theory given above. Interestingly, it has been found that the field from a 
vortex is similar to that from a magnetic monopole located a distance zq = — 1.27A (A is the 



penetration depth) below the superconductor surface ||2^. In this approximation the scalar 
potential from a vortex can be written as 

(27) 



2vr/io + y'^ + {z- zqY ' 
where $0 is the flux quantum. 
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To get some insight into the behaviour of the magnetooptic response, let us neglect 
the absorption and the multiple reflections, and assume that Ha 3> Hr, H^. Then the 
off-diagonal components of the reflection becomes 

Tps oc Qtoi^io [^v{d) -^v{d + D)] + DtanOisincpM^ , (28) 

and 

rsp (X -Qt'A [-^ [<Pv{d) -<pv{d + D)] - Dtan^ism^Af} , (29) 

where d is the distance between the superconductor and the indicator. Although these 
expressions are approximations, they still contain some interesting physics. Note that the 
second term is constant, i.e. independent of the vortex field, since we have assumed that 
Ha ^ Hri Hz- It is also seen that the off- axis reflection coefficients are strongly dependent 
on the penetration depth of the superconductor. That is, when the penetration depth is 
large, the magnetic monopole appears to be located far from the surface of the supercon- 
ductor, whereas when it is small it is located close to the surface. Thus the signal increases 
with decreasing penetration depth. As an example of this. Fig. ^ shows (without the 
constant term) as a function of the penetration depth when d=125 nm (dashed line) and 
d=250 nm (solid line). It is noted that the magnetooptical signal falls off nonlinearly with 
the penetration depth, but that the gradient is less severe when the distance between the 
superconductor and the indicator increases. 

To find the actual signal at the array detector, one must evaluate Eq. |2^. When the 
polarizer and analyzer are in nearly crossed positions (a ~ 90°), and 0c = 0°, the signal can 
be approximated by 

in oc g2 _ ^^(^d + D)f\ r t'tlQsinecoseJo{kr'sine)de\^ . (30) 

L J JO 

Here [(/)^((i) — (t)y^{d + D)]^ may be interpreted as the geometric image of a given position 
(r, 0) at the indicator. According to the Huygens-Fresnel principle, each such point generates 
a spherical secondary wavelet. Since these wavelets must pass through the apertures of the 
optical system, the 'perfect' geometric image is smeared out at the image plane. That 
is, each point in the geometric image must be multiplied by a diffraction integral and the 
point becomes a distribution at the detector. One should also note that in the present 
approximation the diffraction integral is modified by the off-axis reflection coefficient through 
the factor tgi^io- this paper we have not considered the response of the array detector, as 
this is assumed to be ideal. 

It is very important to note that the signal is proportional to the difference in magnetic 
potential at the two surfaces of the indicator. For this reason the signal increases with 
thickness, but the signal-distribution also broadens significantly. This issue will be of par- 
ticular importance when the vortices are placed close together, and the flux-profiles overlap. 
Thus one may expect a thin indicator to provide better signal contrast than a thick one, 
in particular since the contrast is often more important than the strength of the signal. To 
illustrate this point. Fig. ^ shows [0y"((i) — 0v"((i + -D)]^ at the detector plane when the 
vortices are placed 2 yum apart, and the thickness of the indicator is D=5 /im (dashed line) 
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and D=l fim (solid line). Although the signal from the thickest films is 1.7 times stronger 
than that of the thinnest film, we see that the contrast is better in the latter case. 

Finally, it is important to point out that in deriving Eq. (^0]) we have done several 
assumptions in order to obtain a simple expression which is readily interpreted. If a detailed 
analysis is required, one should adopt the full formalism presented in this paper. It is our 
hope that this formalism could be useful for future studies on magnetooptic visualization. 

V. CONCLUSION 

A theoretical model for investigating magnetooptic imaging of weak magnetic fields has 
been presented. It takes into account the details of the indicator and the magnetooptic 
system. To demonstrate the applicability of the model, we have tried to understand some 
issues concerning imaging of single vortices in superconductors. It is found that the useful 
signal is strongly dependent on the penetration depth of the superconductor, and also the 
sensitivity and thickness of the indicator. 
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FIGURES 



FIG. 1. The basic setup for magnetooptical imaging. 

FIG. 2. The magnetization vector in the indicator. 

FIG. 3. The off-axis reflection coefficient as a function of the penetration depth when d=125 
nm (dashed line) and d=250 nm (sohd hne). Here D=l /xm. 

FIG. 4. The normalized signal from three vortices placed 2 /xm apart. Here A = 100 nm, 
D=l ij,m (solid line) and D=5 /xm (dashed line). We have neglected diffraction and assumed a 
magnification M=l. 
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